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We give a simple, semianalytic theory of the development of the hadronic and neu- 

trino/muon components of a cascade induced by a primary which produces hadrons at 

, the initial interaction. The main purpose of the theory is to allow the user to obtain 

Q>^ I quick, but reasonably reliable estimates of the longitudinal properties of such cascades 

developing in a medium, such as a stellar interior, the atmosphere and/or water. As 
an application, we discuss the possibility of discovering physics beyond the Standard 

^^ , Model by means of neutrino telescopes. Some of those events may have spectacular 

^ I signatures in neutrino telescopes. 

:a. 

I ■ 1 Introduction 

D , 
^ . "Back of the envelope" estimates play an important role in physics. First, if 

one wants to decide the feasibility of an experiment, or the observability of a 

K/i ; phenomenon a theorist just discovered, typically, one does not want to spend hours 

}-{ • of CPU time on a substantial computer in order to learn the answer. (It is "no" 

in 90% of the cases, anyway. . . ) Second, even after an elaborate computation has 

been performed, one usually wants to obtain an intuitive understanding of the 

result; in part, in order to learn how to proceed, in part, perhaps, in order to 

discover some bug in the program the presence of which is not obvious otherwise. 

Cascade theory, necessary in order to conduct feasibility studies for neutrino 

telescopes or to understand the generation of neutrinos inside the point sources 

in the sky, among other things, is notoriously time consuming and unintuitive if 

it is to be sufficiently accurate. It is desirable therefore to design some simple 

"back of the envelope" cascade theory for the purpose outlined above. In the age 

before large scale computation became feasible, cascade theory was full of — often 

unreliable — approximations and clumsy numerical computations (cf. the review 

^Invited talk given at the third NESTOR workshop, Pylos 1993; to be published in the 
Proceedings 
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of Nishimura|l]]). This was rendered largely obsolete by the arrival of high speed 
computers. Nevertheless, one can learn some tricks from the "classics" (for good 
reviews, see e.g. Rossi's book, or Nishimura's review article just cited) and 
invent a few new ones in order to reduce the amount of numerical computation 
necessary for obtaining quick estimates. The purpose of the work reported here is 
a first step in this direction; a similar approach to the problem was taken recently 
by Lipari 0. 

2 Approximations 

We outline the approximations used in developing the theory. Most of them can 
be improved upon at some cost in computing time. However, even in its present 
form, the theory is suitable for, say, astrophysical calculations, where uncertainties 
in the input data {e.g. stellar structure, etc. ) are substantially bigger than the 
errors introduced by the approximations. 

• Throughout this work, we use a one dimensional cascade theory in the dif- 
fusion approximation. Transverse development can be added on later in the 
diffusion approximation with relative ease. However, the present version 
is adequate for purposes of a first orientation on ultra high energy (UHE) 
processes. 

• We assume the validity of Feynman scaling in the one particle inclusive cross 

sections, viz. 

da f E \ 

where E and Ein stand for the observed and initial particle energies, respec- 
tively and a stands for the total inelastic cross section. 

In practice, a two parameter fit to F [z] viz. 

F{z) = A{l-zYQ{l~z) (2) 

gives a fair approximation with n ?« 3, see, e.g. |Q. The distribution needs 
an infrared cutoff, in order to get a finite total multiplicity. We choose a 
finite cutoff in z in order to maintain Feynman scaling and thus simplify the 
solution. In order to determine the value of the infrared cutoff, somewhat 
arbitrarily, we chose a median shower energy. El ~ 5 PeV, corresponding 
to a CMS energy, a/s ~ 3 TeV. At this energy, the total average multiplicity 
is (A^) ~ 30, cf.|^. One knows that the infrared cutoff is of the order of 
Aqcd/v^ which, with these numbers gives zq ^ 10~^; in turn, that leads to 
A ^ 4.1. (Most of our results are not very sensitive to the precise value of 
Zq). 



It is known that Feynnian scaling is violated due to QCD loop effects; in 
particular, the average multiplicity of a produced hadron, 

{N) = f-F [z) (3) 

increases with £^j„. However, all such violations of Feynnian scaling have a 
logarithmic dependence on the energy; within the accuracy of the present 
calculations it is legitimate to neglect them. For the sake of consistency, one 
should then either take all logarithmic dependences into account or none of 
them. 

The total inelastic ttN and NN cross sections differ by about 30% or so 
at high energies. Most of the secondaries are pions, with pions of either 
charge being produced in roughly equal numbers. There are few mesons 
produced containing s, c, b (and t ?) quarks. The qualitative reason is 
that those quarks have to be pulled out of a quark sea containing a small 
fraction of heavy quarks. Likewise, baryon pair production is suppressed 
due to the fact that a pair of three quarks has to be created coherently. 
As a consequence, the total high energy baryon contents of a high energy 
interaction is small: most of the baryons (with the exception of a leading 
baryon in a NN interaction) come from target fragmentation and they are of 
low energy in the LAB system. These qualitative arguments are generally 
borne out by the experimental data, wherever they are available, see e.g. f^. 
In view of the above, it is justified to use only one hadron distribution 
function as a first approximation. Instead of a coupled system of integro- 
differential equations, we now have a single equation describing the evolution 
of the cascade. Due to the linearity of the eqations, whatever method of 
solution is used for the single equation, it can be immediately generalized 
to the more accurate, coupled system by a straightforward use of matrix 
methods. However, the increase in computing time is not negligible, since 
one has to invert matrices repeatedly. 

Neutral pions drop out from the hadronic part of the cascade development 
due to their short lifetime at all energies of interest. Therefore, it is justified 
to multiply the hadron distribution by a factor of (2/3) in the cascade evo- 
lution: due to Pomeranchuk's theorem, the vr^N cross sections being very 
nearly equal. Conversely, in the development of the electromagnetic com- 
ponent, the source of photons is fairly represented by 1/3*^ of the hadron 
distribution. (The electromagnetic component of a cascade in this approxi- 
mation will be described elsewhere.) 

Keeping the needs of neutrino telescopes in mind, it is fair to say that a 
neutrino-rich shower is necessarily a hadron-rich one: the main source of 
UHE neutrinos is the process vr^ -^ yU + z/^, where we have not distinguished 



between z/^ and I^^. (The decay of muons is another source of neutrinos; 
however, if one concentrates on the the UHE part of the neutrino spectrum, 
/x-decay neutrinos represent only a small correction, see 0.) Consequently, 
if the initial interaction produces hadrons in substantial numbers, either 
because the primary is a hadron or because of some "new physics" , such as 
the one conjectured in ref. [^ or 't Hooft's B+L violating process, cf. A. 
Ringwald's contribution to these Proceedings, there will be a substantial 
number of neutrinos present. Likewise, if one wants to determine, say, the 
UHE neutrino spectrum emerging from an AGN or a binary system of stars, 
one can concentrate on the accelerated hadrons (p, n. He, etc. ) interacting 
with the target material, which often consists of hadrons. 

In turn, this circumstance drastically simplifies the description of the evo- 
lution of the cascade. Due to the fact that 

(t(7, hadron) ^ Kaa{hadron, hadron), (4) 

(a being the fine structure constant and K a number of 0(1)), the feed- 
back of the electromagnetic component into the hadronic one represents a 
correction of about 1% to the hadronic development. Neutrino-hadron in- 
teractions have a much smaller cross section, typically scaled down by an 

additional factor of ^j^/?7igauge ; where the masses involved are a typical 
hadron mass (say, 1 GeV) and a gauge boson (W, Z) mass, of the order of 
100 GeV. 

Consequently, the evolution of the hadronic component is, in essence, an au- 
tonomous one. For the purposes of determining the neutrino (and electron- 
photon) spectrum, the hadronic component acts only as a source evolving 
autonomously. 

3 The Development of the Cascade 

We denote the differential distribution of hadrons by H{E,x), where x stands 
for the depth measured in units of the hadronic interaction mean free path. In 
this manner, in the UHE region, one can adjust the units as new data become 
available. With current data extrapolated to CMS energies of the order of 50 to 
100 TeV and averaging over vrN and NN cross sections, one gets A ~ 30g/cm^. 
Using the approximations described in the previous Section, we get the diffusion 
equation: 

Here the quantity D stands for the loss of charged hadrons due to decay; its 

expression is: 

Am , , 



where p [x) stands for the density of the medium expressed as a function of a]^. To 
a good approximation, m can be replaced by the pion mass and r by the charged 
pion hfetime, i.e. r ^ 7.8 meters. 

Clearly, the term D is quite small in the UHE region and at a first approach, 
it may be neglected. 

The expression of p (x) depends on the model of the medium used. For in- 
stance, assuming a radial incidence of the primary, we have for some typical 
models: 



1. Exponential atmosphere: 



p (/) = po exp 



/ 



(7) 
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2. Atmosphere with a power law for the density:: 



(8) 



P = Pi 
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p (x) = pi 



x\ 
pihi 



(10) 



3. For a (nearly) incompressible medium, like water, one may even contemplate 
replacing p by a constant, say, the value of the density halfway between the 
entrance of the primary and the depth of the detector. 

For a non-radial incidence, / is to be replaced by the slant depth: the changes in 
the formulae are obvious and will not be exhibited here. 



4 Solutions 



4.1 General techniques 

Neglecting the decay term in eq. (|^), — an approximation valid at the highest 
energies — one can solve the equation by means of a Mellin transformation in a 
standard fashion. The introduction of a new function, h [E, x) by means of the 
substitution 

H {E, x) = exp (— x) h {E, x) 

■^Throughout this work we use natural units, i.e. fi = c = 1. 



reduces the equation to the form: 

dh {E, x) r'^ dE' fE\ 



dx 



I ir^ ye (£'-£) MB'.-). (11) 



This equation is to be solved with an initial condition suitable for the physical 
problem at hand: we outline two types of such problems. 

1. Problems related to the development of a shower in an astrophysical envi- 
ronment or a hadronic shower in the atmosphere. 

In that case, the primary spectrum can be well approximated by a power 
spectrum for a substantial energy range. However, a power spectrum over 
all energies is unphysical, since it has an infinite energy contents^. (We 
also remark in passing that a pure power spectrum has no Mellin transform 
either; thus, the popular factorized solution, (see e.g. 0, p or 0) cannothe 
obtained by using Mellin transforms, just by substituting the ad hoc Ansatz, 
h = E'^'fi^x) into eq. (pA]).) One has to introduce at least an ultraviolet 
cutoff, reflecting the fact that no physical system is capable of producing 
particles of arbitrarily high energy. Often, an infrared cutoff is also needed; 
however, in the present case, this is not necessary: we are concentrating on 
the UHE part of the cascade. Thus, energy-momentum conservation alone 
provides an effective infrared cutoff. 

In the following examples we use a spectrum of the form: 

H{E,0) = n(^) '^QiEM-E) (12) 

2. Problems connected with exploring some "new physics". In all models of 
physics beyond the Standard Model, one has to take into account the fact 
that the Standard Model is highly accurate up to LEP and TEVATRON en- 
ergies. As a consequence, it is commonly assumed (for instance in 0) that 
the onset of the "new physics" is sudden and that it can be characterized 
by a "pseudothreshold" in px and/or energy. Translated to the language 
of cascade development, this means that all but the first interaction is de- 
scribed, in essence, by the Standard Model. Consequently, the effect of the 
"new physics" can be simulated by an initial condition imposed at the depth 
of the first interaction, 

Xq = — — , (13) 



where Xnew is the mfp. corresponding to the new physics introduced. (In 
practice, one may impose the initial condition at depth zero, but x has to 
be replaced by x — Xq.) 



^It is the consequence of this infinite "energy reservoir" that a pure power spectrum passes 
through an absorber without getting distorted. 



There is an important consequence of this picture, often ignored in the lit- 
erature. Some authors have claimed that processes hke 't Hooft's instanton 
induced B+L violating process and other processes, presently considered 
exotic ones, can be recognized by searching for the occurrence of high mul- 
tiplicity muon bundles in neutrino detectors, roughly resembling a heavy 
nucleus hitting the target (the atmosphere or water). While this appears 
to be true for, say, u and Pb induced showers started at the same depth, in 
practice, Xq ^ 1 in any model available: thus, care is needed in counting 
the multiplicity of muon bundles in a neutrino telescope in order to find the 
onset of "new physics" . . . 

With this, the solution of eq. (^ is given by 

H {E, x) = e"^'-^ / ds E~'H (s, 0) exp \f {s - 1) xl , (14) 

where H and F stand for the Mellin transforms of the initial condition and of F, 
respectively. The contour C runs from —zoo to -\-ioo in the complex s plane, in 
the strip where both F and H exist. 

Equation (|T^) can be a starting point of a numerical evaluation of the hadronic 
distribution. The initial conditions mentioned above (and several variations on 
them) as well as the expression of F are simple enough so that their Mellin trans- 
forms can be computed analytically; yet, they are sufficiently accurate so as to 
give a fair representation of the initial spectra and of the inclusive cross sections. 
In general, however the integral in eq. ([T^) cannot be evaluated in a closed form. 
Asymptotic methods used in the distant past (saddle point, etc. ) are, in general, 
not sufficiently accurate. 

Another way, very convenient from the point of view of a numerical treatment, 
is to solve the diffusion equation by successive approximations. The convenient 
starting point for this is eq. ([TT|). Upon integrating both sides with respect to x 
and putting the equation in the form of a recursion relation, we obtain: 

h^^+'^ {E, x) = h {E, 0) + J^ j^ dx'—F (^— j /i(") {E', x') . (15) 

Here /i*-"^ stands for the n*'^ iteration of the solution, with 

h^ = h{E,0). 

One notices that the integration over the variable x' can be performed in a 
closed form at every step of the iterative procedure. Moreover, with the simple 
representation of the inclusive distribution suggested in this section and with the 
simple initial spectra descibed above, the integral over the energy can also be 
evaluated in a closed form, although the result is somewhat clumsy. In this way, 
the iterative solution reduces to the evaluation of a series. ( This, however, may 
not be true for more complicated representations of the initial inclusive distribution 
and/or of the initial spectrum.) 

A few remarks are in order here. 



1. One knows that the iterative process just described is a convergent one: 
eq. ([T5|) is a Volterra equation which (due to the compactness of its kernel) 
can always be solved by successive approximations. (For purposes of prov- 
ing this fact, one should use 1/E' as an independent variable in eq. (p^)). 
As a consequence, if one uses the simple representations of the inclusive 
cross section and initial spectra as described above, the resulting series is a 
convergent one (in practice, it converges quite rapidly). 

2. The method of successive approximations described above is, in essence, 
equivalent to the method of successive collisions invented by Bhabha and 
Heitler in 1937, see e.g. [0. 

3. It is often advantageous to combine the numerical evaluation of the cascade 
development with symbolic manipulation programs, such as MATHEMAT- 
ICA(C). In our experience, often a substantial amount of computing time 
may be saved in this manner. 

4.2 A Very Simple Solution 

If one is interested in the crudest type of estimates only, one does not worry too 
much about the precise shape of the inclusive distribution. In particular, one may 
just assume that there is no leading particle present at all and the energy is shared 
equally among all the secondaries, i.e. 
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1 



(N), 



(16) 



The validity of this approximation has been discussed elsewhere, see . We 
concluded in that reference that the relative error committed by using eq. ([I6| ) as 
opposed to a more realistic form of the inclusive distribution is about a factor of 
2. Hence, even this extremely crude approximation is adequate for the purposes 
of qualitative estimates. Using eq. (|T6|) , one easily arrives at an explicit solution. 
Assuming a primary power spectrum cut off at i? = Em, one gets: 



H{E,x) =E-''e 



k=0 



2-a\^ 



k! 



(17) 



The maximal value of k is given by: 



E 



M 



E\n{N) 



Here [■ ■ ■] stands for the integer part of a number. 

Due to the sharp cutoff of the spectrum, H {E, x) is a function which has 
discontinuities. On could get rid of those discontinuities by smoothing out in some 
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way the B function occurring in the spectrum: in that case, however, one would 
be deahng with an infinite series. Due to the crudeness of this approximation, 
we shall not discuss this question any further. One notices that the approximate 
solution written down here is a simple generalization of "Heitler's caricature" of 
an electromagnetic cascade, cf. 0. 

4.3 Muons and Muon Neutrinos 

In the approximation used in the present work (UHE neutrinos and muons only, no 
decaying muons, cf. the previous Section), there are practically no u^, present: the 
branching ratio of vr^ — >■ e + 1/ is of the order of 10~^, cf. 0. Further, the number 
of neutrinos and charged muons is equal and determined by the autonomously 
evolving hadronic component: 

N{u + V)=N(fi+ + fi-) =^ N. 



We get in an obvious manner: 

N {E, x) = - rdx' D {2E, x') H {2E, x') 
3 JO 



Thus, once the hadronic component is computed by using one of the solu- 
tions described in the previous subsections, the number of UHE muons and muon 
neutrinos is obtained by a quadrature. (It is known of course that, either in the 
case of a pure power spectrum, or going to lower energies in the present calcula- 
tion, electron neutrinos begin to appear. At energies substantially lower than the 
ultraviolet cutoff, the ratio of muon and electron neutrinos is about 2.) 

Given the simplifications made in the course of the calculations, one should 
check the accuracy and internal consistency of the results. 

First, an internal consistency check. We stated earlier that the iterative solu- 
tion to the cascade equation should be convergent on general grounds. In practice, 
however, the speed of the convergence is an important factor from the computa- 
tional point of view. For this reason, we computed the hadronic fiux and the 
atmospheric neutrino flux generated by it at a depth of 1000 g/cm^ and with a 
primary spectrum oc i?~" cut off at Ej^j = lO^^GeV. For the sake of definiteness, 
we chose a = 2.7. 

The iterative solution can be written in the form: 

oo 

H{x, E) oc E-^e-^ ^ (ax)" /„ {E/Eu) , (19) 

n=0 

where the constant a is determined in terms of fractional moments of the frag- 
mentation function. One easily verifies that as Em — >■ cc, all /„ — *> 1 and the 
solution goes over into the one obtained with the factorization assumption. 

In order to exhibit the speed of the convergence, we chose /o = 1. In Fig. |I] 
we exhibited the first few coefficients fn{E / Em)- 
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Figure 1: The first six coefficients of tlie iterative solution. Tlie consecutive 
coefficients are plotted on a varying gray scale, /i being the darkest and widest 
line. 

One sees that the iterative solution is converging quite well: the area under 
the consecutive coefficients fn{E j Em) decreases rapidly. 

Next, we compute the hadron distribution and from it, the neutrino flux in the 
manner described above. For the sake of simplicity, we used a standard expression 
of the flux, ref. p] all the way to the cutoff. In Fig ^ we plot the integral spectrum 
of neutrinos, multiplied by E'^''^ . 



If the factorized solution were valid everywhere, this quantity would be ap- 
proximately constant in the energy range we are considering. 

We see that this is reasonably well satisfied up to E/Em ~ 10^^. In this 
region, our result agrees to within 30% or so, with Lipari's0, and also by MC 
calculations. Given the simplicity of the model developed here, the agreement is 
satisfactory. 

At still higher energies, the effect of the ultraviolet cutoff can be clearly seen. In 
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Figure 2: The integral spectrum of atmospheric neutrinos at a depth 1000 g/cm^. 

principle, one could infer the existence of such a cutoff from atmospheric neutrino 
observations. The main obstacle to such an observation is, of course, the scarceness 
of events. 

5 Looking for "New Physics" 

In this Section, we illustrate the method developed in the previous Sections by 
asking a simple question: 

Assuming that there exists some "new physics" beyond the Standard Model, do 
we have any hope of detecting it in a neutrino telescope ? 

The answer, of course, depends on the cross section and energy of the onset 
of the new processes. We have in mind, in particular, either the phenomenon 
conjectured in ref.0 or the multiple production of gauge Bosons (with or without 
B+L violation) as discussed by A. Ringwald in these Proceedings. 

There is a large amount of theoretical uncertainty in the nature of the phe- 
nomena in which any physics beyond the Standard Model would manifest itself. 
However, one can determine with relative ease the reactions in which it is virtually 
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hopeless to look for manifestations of some post-Standard-Model physics, unless 
the effects are unexpectedly dramatic. Any such argument is based, in essence, 
on the unitarity of the S-matrix. Probability is conserved, hence, one should not 
look for the new physics in reactions with a large number of open channels. Unless 
a firm theoretical prediction is available — as it was in the case of the discovery 
of the W and Z bosons in a hadronic machine — it is virtually impossible to find 
the proverbial needle in a haystack. . . 

This criterion immediately tells us that hadronic reactions are, in all probabil- 
ity, unsuitable in a search for new physics: there are just too many channels open, 
producing mostly mundane physics: typically, relatively soft pions. Likewise, 
electron pair production in 7 - nucleus collisions dominates in a medium with an 
average Z of the order of 3 or larger (water, earth, the atmosphere). Hence, again, 
"mundane physics" suppresses any appearence of new physics. (This question 



was discussed in some detail in ref. ||T0[- In a recent work, Morris and Ringwald, 
ref. [^ reached a similar conclusion by means of a rather sophisticated Monte 
Carlo simulation.) 

What remains therefore, is the realm of neutrino induced reactions: any pro- 
cess predicted by the Standard Model has a cross section of the order of a few 
nanobarns at LAB. energies around an EeV or so. By contrast, the conjectured 
new processes may reach cross sections of the order of 10~^mb, cf. ref. 0, [p!T[| . 
Therefore, they should be observable in neutrino induced reactions better than 
anywhere else. 

Neglecting all effects of nuclear structure, there is a very simple relationship 
between a mfp (A) and the cross section of a reaction {cf. and references 
discussed there): 

Afe/cT ^ i^. (20, 

(At high energies, this relationship should hold to an accuracy better than about 
20% for medium heavy nuclei.) 

Thus, a //-induced reaction of a cross section, cr ^ 5 x 10~'^mb and incidence 
at a zenith angle, 6 = 0, would, in the mean, produce its first interaction very 
close to a detector like NESTOR. (At incidences lower than vertical, the range 
of cross sections one can explore below the one at vertical incidence, depends on 
the geology around the detector and it needs a more detailed investigation. In 
general, however, one can explore several orders of magnitude in the initial cross 
section by scanning at lower zenith angles. For purposes of illustration therefore, 
we assume vertical incidence and a cross section of the order of magnitude just 
quoted.) 

Due to the theoretical uncertainties surrounding the "new physics" , we made 
two simplifying assumptions. 

1. The onset of the "new physics" is sudden: in practice, we approximate it by 
a function in the energy, as in the articles in ref. [|^. As a consequence, 
after the initial interaction, the resulting shower tends to evolve according 
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to the Standard Model. (At extremely high energies where this is not the 
case, the primary fluxes are expected to be very low.) 

2. The rapidity distribution of the produced particles in the initial interaction 
follows an "equipartition law", viz. 

F(z) = ^ (^ - ^) . (21) 

where A^^i stands for the average hadronic multiplicity in the first interaction: 
this quantity parametrizes what we call "new physics" in this model. 

Rather than compute the spectrum of the hadronic component for a given 
primary neutrino spectrum, we investigated the profile of single events, with a 
primary energy E = 10^^ GeV, corresponding to about ^/s ~ 5 TeV in the 
neutrino-nucleon CMS. 

We computed the development of the hadronic cascade with two initial multi- 
plicities (A^i = 5 and A^i = 20). 

Both cascades are supposed to have a primary incident at zenith angle, 9 = 
and an initial cross section around the value just quoted. The results are displayed 
in Fig. I and Fig. ||, respectively. 



In both Figures, we plotted the integral hadron spectrum as a function of 
X = t/X for energies E > 10 GeV and E > 100 GeV, respectively. (Obviously, the 
higher curves correspond to the lower threshold energies in both Figures.) 

In the Figures, A is the usual hadronic mfp, given the assumption that after 
the first interaction, the cascade develops according to the Standard Model. 



The interesting feature of these showers is that charged hadrons containing 
light quarks (u,d,s), for all practical purposes, do not decay. (Indeed, the inter- 
action mfp is about 80 cm, to be compared with decay mfps substantially larger 
than 10 meters.) Similarly, the radiation length is approximately 36 cm. Hence, 
both the hadronic and electromagnetic cascades evolve within a few meters in real 
space. The 'anomalous' showers are characterized by a very large energy deposi- 
tion within the span of a few meters. In particular, a large amount of Cherenkov 
light is emitted in such events. 

Some features are worth noticing. 

• The overall profile of the longitudinal evolution of a shower is rather insensi- 
tive to the initial multiplicity: hence, theoretical uncertainties are unlikely to 
have a drastic effect on the discovery potential of a given neutrino telescope. 
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Figure 3: The profile the hadronic component of an 'anomalous' shower with 

iVi = 5. 

• The events are "muon poor" for the reason stated above. In practice, the 
only source of muons is the production of particles containing heavy quarks 
(c, b, t). However, the production of heavy quarks is expected to be sup- 
pressed due to the scarcity of the latter in the quark sea. (In the case of 
multiple W-production, the source of muons is the leptonic decay mode of 
the gauge bosons, which is of the order of 30%.) 

• Overall, one concludes that an underwater neutrino telescope may be the 
detector of choice for the observation of the events discussed here. Due to 
the large local energy deposition, such anomalous events may be observable 
even somewhat beyond the usual attenuation length of Cherenkov light (of 
the order of 50m at NESTOR). This question needs further study. 

The potential beauty of a search for new physics in observations of the type 
just described is that the signature is a very robust one. Events in which a large 
amount of energy is released near the neutrino telescope are hard to miss, even 
though one does not understand all the details of such events at present. 
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Figure 4: The profile of the hadronic component of an 'anomalous' shower with 
A^i = 20. 

One event of an apparently large energy release and multiplicity has been 
reported recently by the KAMIOKANDE collaboration, ref. |]T2|]. It is amusing 



to speculate about the possibility that the "unusual event" as described by the 
authors may be the manifestation of a phenomenon just discussed. Should that 
be the case, one would have the first experimental evidence for the existence of 
some physics by means of which the ills of the Standard Model could be cured. 

6 Discussion 

We believe that the formalism presented in this paper is useful in order to ex- 
plore new phenomena both in an astrophysical context (point sources, AGN. . . ) 
and in a particle physics one. Neutrino telescopes are potentially sensitive to 
the observation oi u + ^quark' reactions at energies beyond the reach of current 
accelerator-based research. Hence, despite the erratic nature of point sources and 
the uncertainties associated with sources like AGN or gamma - bursters, they may 
have an important role to play in exploring particle physics beyond the ^/s ~ 1 
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TeV barrier. 
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